CRESSP: a comprehensive pipeline for prediction of immunopathogenic SARS-CoV-2 epitopes using structural properties of proteins.

The development of autoimmune diseases following SARS-CoV-2 infection, including multisystem inflammatory syndrome, has been reported, and several mechanisms have been suggested, including molecular mimicry. We developed a scalable, comparative immunoinformatics pipeline called cross-reactive-epitope-search-using-structural-properties-of-proteins (CRESSP) to identify cross-reactive epitopes between a collection of SARS-CoV-2 proteomes and the human proteome using the structural properties of the proteins. Overall, by searching 4 911 245 proteins from 196 352 SARS-CoV-2 genomes, we identified 133 and 648 human proteins harboring potential cross-reactive B-cell and CD8+ T-cell epitopes, respectively. To demonstrate the robustness of our pipeline, we predicted the cross-reactive epitopes of coronavirus spike proteins, which were recognized by known cross-neutralizing antibodies. Using single-cell expression data, we identified PARP14 as a potential target of intermolecular epitope spreading between the virus and human proteins. Finally, we developed a web application (https://ahs2202.github.io/3M/) to interactively visualize our results. We also made our pipeline available as an open-source CRESSP package (https://pypi.org/project/cressp/), which can analyze any two proteomes of interest to identify potentially cross-reactive epitopes between the proteomes. Overall, our immunoinformatic resources provide a foundation for the investigation of molecular mimicry in the pathogenesis of autoimmune and chronic inflammatory diseases following COVID-19.

Predicted HLA Class I and Class II Epitopes From Licensed Vaccines Are Largely Conserved in New SARS-CoV-2 Omicron Variant of Concern.

The potential effect of emerging SARS-CoV-2 variants on vaccine efficacy is an issue of critical importance. In this study, the possible impact of mutations that facilitate virus escape from the cytotoxic and the helper cellular immune responses in the new SARS-CoV-2 Omicron variant of concern was analyzed for the 551 and 41 most abundant HLA class I and II alleles, respectively. Computational prediction showed that almost all of these 592 alleles, which cover >90% of the human population, contain enough epitopes without escape mutations in the emerging SARS-CoV-2 Omicron variant of concern. These data suggest that both cytotoxic and helper cellular immune protection elicited by currently licensed vaccines are virtually unaffected by the highly contagious SARS-CoV-2 Omicron variant of concern.

A class II MHC-targeted vaccine elicits immunity against SARS-CoV-2 and its variants.

The pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in over 100 million infections and millions of deaths. Effective vaccines remain the best hope of curtailing SARS-CoV-2 transmission, morbidity, and mortality. The vaccines in current use require cold storage and sophisticated manufacturing capacity, which complicates their distribution, especially in less developed countries. We report the development of a candidate SARS-CoV-2 vaccine that is purely protein based and directly targets antigen-presenting cells. It consists of the SARS-CoV-2 Spike receptor-binding domain (SpikeRBD) fused to an alpaca-derived nanobody that recognizes class II major histocompatibility complex antigens (VHHMHCII). This vaccine elicits robust humoral and cellular immunity against SARS-CoV-2 and its variants. Both young and aged mice immunized with two doses of VHHMHCII-SpikeRBD elicit high-titer binding and neutralizing antibodies. Immunization also induces strong cellular immunity, including a robust CD8 T cell response. VHHMHCII-SpikeRBD is stable for at least 7 d at room temperature and can be lyophilized without loss of efficacy.

Can molecular mimicry explain the cytokine storm of SARS-CoV-2?: An in silico approach.

PARP14 and PARP9 play a key role in macrophage immune regulation. SARS-CoV-2 is an emerging viral disease that triggers hyper-inflammation known as a cytokine storm. In this study, using in silico tools, we hypothesize about the immunological phenomena of molecular mimicry between SARS-CoV-2 Nsp3 and the human PARP14 and PARP9. The results showed an epitope of SARS-CoV-2 Nsp3 protein that contains consensus sequences for both human PARP14 and PARP9 that are antigens for MHC Classes 1 and 2, which can potentially induce an immune response against human PARP14 and PARP9; while its depletion causes a hyper-inflammatory state in SARS-CoV-2 patients.

A Systematic, Unbiased Mapping of CD8+ and CD4+ T Cell Epitopes in Yellow Fever Vaccinees.

Examining CD8+ and CD4+ T cell responses after primary Yellow Fever vaccination in a cohort of 210 volunteers, we have identified and tetramer-validated 92 CD8+ and 50 CD4+ T cell epitopes, many inducing strong and prevalent (i.e., immunodominant) T cell responses. Restricted by 40 and 14 HLA-class I and II allotypes, respectively, these responses have wide population coverage and might be of considerable academic, diagnostic and therapeutic interest. The broad coverage of epitopes and HLA overcame the otherwise confounding effects of HLA diversity and non-HLA background providing the first evidence of T cell immunodomination in humans. Also, double-staining of CD4+ T cells with tetramers representing the same HLA-binding core, albeit with different flanking regions, demonstrated an extensive diversification of the specificities of many CD4+ T cell responses. We suggest that this could reduce the risk of pathogen escape, and that multi-tetramer staining is required to reveal the true magnitude and diversity of CD4+ T cell responses. Our T cell epitope discovery approach uses a combination of (1) overlapping peptides representing the entire Yellow Fever virus proteome to search for peptides containing CD4+ and/or CD8+ T cell epitopes, (2) predictors of peptide-HLA binding to suggest epitopes and their restricting HLA allotypes, (3) generation of peptide-HLA tetramers to identify T cell epitopes, and (4) analysis of ex vivo T cell responses to validate the same. This approach is systematic, exhaustive, and can be done in any individual of any HLA haplotype. It is all-inclusive in the sense that it includes all protein antigens and peptide epitopes, and encompasses both CD4+ and CD8+ T cell epitopes. It is efficient and, importantly, reduces the false discovery rate. The unbiased nature of the T cell epitope discovery approach presented here should support the refinement of future peptide-HLA class I and II predictors and tetramer technologies, which eventually should cover all HLA class I and II isotypes. We believe that future investigations of emerging pathogens (e.g., SARS-CoV-2) should include population-wide T cell epitope discovery using blood samples from patients, convalescents and/or long-term survivors, who might all hold important information on T cell epitopes and responses.

Individual HLA-A, -B, -C, and -DRB1 Genotypes Are No Major Factors Which Determine COVID-19 Severity.

HLA molecules are key restrictive elements to present intracellular antigens at the crossroads of an effective T-cell response against SARS-CoV-2. To determine the impact of the HLA genotype on the severity of SARS-CoV-2 courses, we investigated data from 6,919 infected individuals. HLA-A, -B, and -DRB1 allotypes grouped into HLA supertypes by functional or predicted structural similarities of the peptide-binding grooves did not predict COVID-19 severity. Further, we did not observe a heterozygote advantage or a benefit from HLA diplotypes with more divergent physicochemical peptide-binding properties. Finally, numbers of in silico predicted viral T-cell epitopes did not correlate with the severity of SARS-CoV-2 infections. These findings suggest that the HLA genotype is no major factor determining COVID-19 severity. Moreover, our data suggest that the spike glycoprotein alone may allow for abundant T-cell epitopes to mount robust T-cell responses not limited by the HLA genotype.

Influence of HLA Class II Polymorphism on Predicted Cellular Immunity Against SARS-CoV-2 at the Population and Individual Level.

Development of adaptive immunity after COVID-19 and after vaccination against SARS-CoV-2 is predicated on recognition of viral peptides, presented on HLA class II molecules, by CD4+ T-cells. We capitalised on extensive high-resolution HLA data on twenty five human race/ethnic populations to investigate the role of HLA polymorphism on SARS-CoV-2 immunogenicity at the population and individual level. Within populations, we identify wide inter-individual variability in predicted peptide presentation from structural, non-structural and accessory SARS-CoV-2 proteins, according to individual HLA genotype. However, we find similar potential for anti-SARS-CoV-2 cellular immunity at the population level suggesting that HLA polymorphism is unlikely to account for observed disparities in clinical outcomes after COVID-19 among different race/ethnic groups. Our findings provide important insight on the potential role of HLA polymorphism on development of protective immunity after SARS-CoV-2 infection and after vaccination and a firm basis for further experimental studies in this field.

Identification of presented SARS-CoV-2 HLA class I and HLA class II peptides using HLA peptidomics.

The human leukocyte antigen (HLA)-bound viral antigens serve as an immunological signature that can be selectively recognized by T cells. As viruses evolve by acquiring mutations, it is essential to identify a range of presented viral antigens. Using HLA peptidomics, we are able to identify severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-derived peptides presented by highly prevalent HLA class I (HLA-I) molecules by using infected cells as well as overexpression of SARS-CoV-2 genes. We find 26 HLA-I peptides and 36 HLA class II (HLA-II) peptides. Among the identified peptides, some are shared between different cells and some are derived from out-of-frame open reading frames (ORFs). Seven of these peptides were previously shown to be immunogenic, and we identify two additional immunoreactive peptides by using HLA multimer staining. These results may aid the development of the next generation of SARS-CoV-2 vaccines based on presented viral-specific antigens that span several of the viral genes.

Mapping the SARS-CoV-2 spike glycoprotein-derived peptidome presented by HLA class II on dendritic cells.

Understanding and eliciting protective immune responses to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an urgent priority. To facilitate these objectives, we profile the repertoire of human leukocyte antigen class II (HLA-II)-bound peptides presented by HLA-DR diverse monocyte-derived dendritic cells pulsed with SARS-CoV-2 spike (S) protein. We identify 209 unique HLA-II-bound peptide sequences, many forming nested sets, which map to sites throughout S including glycosylated regions. Comparison of the glycosylation profile of the S protein to that of the HLA-II-bound S peptides reveals substantial trimming of glycan residues on the latter, likely induced during antigen processing. Our data also highlight the receptor-binding motif in S1 as a HLA-DR-binding peptide-rich region and identify S2-derived peptides with potential for targeting by cross-protective vaccine-elicited responses. Results from this study will aid analysis of CD4+ T cell responses in infected individuals and vaccine recipients and have application in next-generation vaccine design.

Identification of SARS-CoV-2 CTL epitopes for development of a multivalent subunit vaccine for COVID-19.

An immunoinformatics-based approach was used to identify potential multivalent subunit CTL vaccine candidates for SARS-CoV-2. Criteria for computational screening included antigen processing, antigenicity, allergenicity, and toxicity. A total of 2604 epitopes were found to be strong binders to MHC class I molecules when analyzed using IEDB tools. Further testing for antigen processing yielded 826 peptides of which 451 were 9-mers that were analyzed for potential antigenicity. Antigenic properties were predicted for 102 of the 451 peptides. Further assessment for potential allergenicity and toxicity narrowed the number of candidate CTL epitopes to 50 peptide sequences, 45 of which were present in all strains of SARS-CoV-2 that were tested. The predicted CTL epitopes were then tested to eliminate those with MHC class II binding potential, a property that could induce hyperinflammatory responses mediated by TH2 cells in immunized hosts. Eighteen of the 50 epitopes did not show class II binding potential. To our knowledge this is the first comprehensive analysis on the proteome of SARS-CoV-2 for prediction of CTL epitopes lacking binding properties that could stimulate unwanted TH2 responses. Future studies will be needed to assess these epitopes as multivalent subunit vaccine candidates which stimulate protective CTL responses against SARS-COV-2.

Assessment of SARS-CoV-2 specific CD4(+) and CD8 (+) T cell responses using MHC class I and II tetramers.

The success of SARS-CoV-2 (CoV-2) vaccines is measured by their ability to mount immune memory responses that are long-lasting. To achieve this goal, it is important to identify surrogates of immune protection, namely, CoV-2 MHC Class I and II immunodominant pieces/epitopes and methodologies to measure them. Here, we present results of flow cytometry-based MHC Class I and II QuickSwitchTM platforms for assessing SARS-CoV-2 peptide binding affinities to various human alleles as well as the H-2 Kb mouse allele. Multiple SARS-CoV-2 potential MHC binders were screened and validated by QuickSwitch testing. The screen included 31 MHC Class I and 19 MHC Class II peptides predicted to be good binders by the IEDB web resource provided by NIAID. While several predicted peptides with acceptable theoretical Kd showed poor MHC occupancies, fourteen MHC class II and three MHC class I peptides showed promiscuity in that they bind to multiple MHC molecule types. In addition to providing important data towards the study of the SARS-CoV-2 virus and its presented antigenic epitopes, the peptides identified in this study can be used in the QuickSwitch platform to generate MHC tetramers. With those tetramers, scientists can assess CD4 + and CD8 + immune responses to these different MHC/peptide complexes.

Immuno-informatics design of a multimeric epitope peptide based vaccine targeting SARS-CoV-2 spike glycoprotein.

Developing an efficacious vaccine for SARS-CoV-2 infection is critical to stemming COVID-19 fatalities and providing the global community with immune protection. We have used a bioinformatic approach to aid in designing an epitope peptide-based vaccine against the spike protein of the virus. Five antigenic B cell epitopes with viable antigenicity and a total of 27 discontinuous B cell epitopes were mapped out structurally in the spike protein for antibody recognition. We identified eight CD8+ T cell 9-mers and 12 CD4+ T cell 14-15-mer as promising candidate epitopes putatively restricted by a large number of MHC I and II alleles, respectively. We used this information to construct an in silico chimeric peptide vaccine whose translational rate was highly expressed when cloned in pET28a (+) vector. With our In silico test, the vaccine construct was predicted to elicit high antigenicity and cell-mediated immunity when given as a homologous prime-boost, triggering of toll-like receptor 5 by the adjuvant linker. The vaccine was also characterized by an increase in IgM and IgG and an array of Th1 and Th2 cytokines. Upon in silico challenge with SARS-CoV-2, there was a decrease in antigen levels using our immune simulations. We, therefore, propose that potential vaccine designs consider this approach.

Population-Predicted MHC Class II Epitope Presentation of SARS-CoV-2 Structural Proteins Correlates to the Case Fatality Rates of COVID-19 in Different Countries.

We observed substantial differences in predicted Major Histocompatibility Complex II (MHCII) epitope presentation of SARS-CoV-2 proteins for different populations but only minor differences in predicted MHCI epitope presentation. A comparison of this predicted epitope MHC-coverage revealed for the early phase of infection spread (till day 15 after reaching 128 observed infection cases) highly significant negative correlations with the case fatality rate. Specifically, this was observed in different populations for MHC class II presentation of the viral spike protein (p-value: 0.0733 for linear regression), the envelope protein (p-value: 0.023), and the membrane protein (p-value: 0.00053), indicating that the high case fatality rates of COVID-19 observed in some countries seem to be related with poor MHC class II presentation and hence weak adaptive immune response against these viral envelope proteins. Our results highlight the general importance of the SARS-CoV-2 structural proteins in immunological control in early infection spread looking at a global census in various countries and taking case fatality rate into account. Other factors such as health system and control measures become more important after the early spread. Our study should encourage further studies on MHCII alleles as potential risk factors in COVID-19 including assessment of local populations and specific allele distributions.

In silico analysis suggests less effective MHC-II presentation of SARS-CoV-2 RBM peptides: Implication for neutralizing antibody responses.

SARS-CoV-2 antibodies develop within two weeks of infection, but wane relatively rapidly post-infection, raising concerns about whether antibody responses will provide protection upon re-exposure. Here we revisit T-B cooperation as a prerequisite for effective and durable neutralizing antibody responses centered on a mutationally constrained RBM B cell epitope. T-B cooperation requires co-processing of B and T cell epitopes by the same B cell and is subject to MHC-II restriction. We evaluated MHC-II constraints relevant to the neutralizing antibody response to a mutationally-constrained B cell epitope in the receptor binding motif (RBM) of the spike protein. Examining common MHC-II alleles, we found that peptides surrounding this key B cell epitope are predicted to bind poorly, suggesting a lack MHC-II support in T-B cooperation, impacting generation of high-potency neutralizing antibodies in the general population. Additionally, we found that multiple microbial peptides had potential for RBM cross-reactivity, supporting previous exposures as a possible source of T cell memory.

T and B cell Epitope analysis of SARS-CoV-2 S protein based on immunoinformatics and experimental research.

COVID-19 caused by SARS-CoV-2 is pandemic with a severe morbidity and mortality rate across the world. Despite the race for effective vaccine and drug against further expansion and fatality rate of this novel coronavirus, there is still lack of effective antiviral therapy. To this effect, we deemed it necessary to identify potential B and T cell epitopes from the envelope S protein. This can be used as potential targets to develop anti-SARS-CoV-2 vaccine preparations. In this study, we used immunoinformatics to identify conservative B and T cell epitopes for S proteins of SARS-CoV-2, which might play roles in the initiation of SARS-CoV-2 infection. We identified the B cell and T cell peptide epitopes of S protein and their antigenicity, as well as the interaction between the peptide epitopes and human leucocyte antigen (HLA). Among the B cell epitopes, 'EILDITPCSFGGVS' has the highest score of antigenicity and great immunogenicity. In T cell epitopes, MHC-I peptide 'KIADYNYKL' and MHC-II peptide 'LEILDITPC' were identified as high antigens. Besides, docking analysis showed that the predicted peptide 'KIADYNYKL' was closely bound to the HLA-A*0201. The results of molecular dynamics simulation through GROMACS software showed that 'HLA-A*0201~peptide' complex was very stable. And the peptide we selected could induce the T cell response similar to that of SARS-CoV-2 infection. Moreover, the predicted peptides were highly conserved in different isolates from different countries. The antigenic epitopes presumed in this study were effective new vaccine targets to prevent SARS-CoV-2 infection.

The Human Leukocyte Antigen Class II Immunopeptidome of the SARS-CoV-2 Spike Glycoprotein.

Precise elucidation of the antigen sequences for T cell immunosurveillance greatly enhances our ability to understand and modulate humoral responses to viral infection or active immunization. Mass spectrometry is used to identify 526 unique sequences from the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike glycoprotein extracellular domain in a complex with human leukocyte antigen class II molecules on antigen-presenting cells from a panel of healthy donors selected to represent a majority of allele usage from this highly polymorphic molecule. The identified sequences span the entire spike protein, and several sequences are isolated from a majority of the sampled donors, indicating promiscuous binding. Importantly, many peptides derived from the receptor binding domain used for cell entry are identified. This work represents a precise and comprehensive immunopeptidomic investigation with the SARS-CoV-2 spike glycoprotein and allows detailed analysis of features that may aid vaccine development to end the current coronavirus disease 2019 (COVID-19) pandemic.

Exploring the out of sight antigens of SARS-CoV-2 to design a candidate multi-epitope vaccine by utilizing immunoinformatics approaches.

SARS-CoV-2 causes a severe respiratory disease called COVID-19. Currently, global health is facing its devastating outbreak. However, there is no vaccine available against this virus up to now. In this study, a novel multi-epitope vaccine against SARS-CoV-2 was designed to provoke both innate and adaptive immune responses. The immunodominant regions of six non-structural proteins (nsp7, nsp8, nsp9, nsp10, nsp12 and nsp14) of SARS-CoV-2 were selected by multiple immunoinformatic tools to provoke T cell immune response. Also, immunodominant fragment of the functional region of SARS-CoV-2 spike (400-510 residues) protein was selected for inducing neutralizing antibodies production. The selected regions' sequences were connected to each other by furin-sensitive linker (RVRR). Moreover, the functional region of ß-defensin as a well-known agonist for the TLR-4/MD complex was added at the N-terminus of the vaccine using (EAAAK)3 linker. Also, a CD4 + T-helper epitope, PADRE, was used at the C-terminal of the vaccine by GPGPG and A(EAAAK)2A linkers to form the final vaccine construct. The physicochemical properties, allergenicity, antigenicity, functionality and population coverage of the final vaccine construct were analyzed. The final vaccine construct was an immunogenic, non-allergen and unfunctional protein which contained multiple CD8 + and CD4 + overlapping epitopes, IFN-γ inducing epitopes, linear and conformational B cell epitopes. It could form stable and significant interactions with TLR-4/MD according to molecular docking and dynamics simulations. Global population coverage of the vaccine for HLA-I and II were estimated 96.2% and 97.1%, respectively. At last, the final vaccine construct was reverse translated to design the DNA vaccine. Although the designed vaccine exhibited high efficacy in silico, further experimental validation is necessary.

COVID-19 in a lung transplant recipient: Exploring the diagnostic role of circulating exosomes and the clinical impact of advanced immunosuppression.

Exosomes isolated from plasma of lung transplant recipients with allograft injury contain donor-derived lung self-antigens (collagen V and Kα1 tubulin) and human leukocyte antigen (HLA) molecules. We present a case of a 76-year-old, female lung transplant recipient treated for acute cellular rejection with methylprednisolone and anti-thymocyte globulin, who subsequently contracted SARS-CoV-2 and developed a sharp increase in the mean fluorescent intensity of anti-HLA antibodies. Analysis of circulating exosomes during rejection, but before SARS-CoV-2 infection, revealed the presence of lung self-antigens and HLA class II molecules. After the patient contracted SARS-CoV-2, exosomes with the SARS-CoV-2 spike protein were also found. After resolution of infectious symptoms, exosomes with SARS-CoV-2 spike protein were no longer detected; however, exosomes with lung self-antigens and HLA class II molecules persisted, which coincided with a progressive decline in spirometric flows, suggesting chronic lung allograft dysfunction. We propose that the analysis of circulating exosomes may be used to detect allograft injury mediated by both rejection and infection. Furthermore, the detection of exosomes containing viral proteins may be helpful in identifying allograft injury driven by viral pathogens.

Computationally Optimized SARS-CoV-2 MHC Class I and II Vaccine Formulations Predicted to Target Human Haplotype Distributions.

We present a combinatorial machine learning method to evaluate and optimize peptide vaccine formulations for SARS-CoV-2. Our approach optimizes the presentation likelihood of a diverse set of vaccine peptides conditioned on a target human-population HLA haplotype distribution and expected epitope drift. Our proposed SARS-CoV-2 MHC class I vaccine formulations provide 93.21% predicted population coverage with at least five vaccine peptide-HLA average hits per person (≥ 1 peptide: 99.91%) with all vaccine peptides perfectly conserved across 4,690 geographically sampled SARS-CoV-2 genomes. Our proposed MHC class II vaccine formulations provide 97.21% predicted coverage with at least five vaccine peptide-HLA average hits per person with all peptides having an observed mutation probability of ≤ 0.001. We provide an open-source implementation of our design methods (OptiVax), vaccine evaluation tool (EvalVax), as well as the data used in our design efforts here: https://github.com/gifford-lab/optivax.

Mortality in COVID-19 disease patients: Correlating the association of major histocompatibility complex (MHC) with severe acute respiratory syndrome 2 (SARS-CoV-2) variants.

Genetic factors such as the HLA type of patients may play a role in regard to disease severity and clinical outcome of patients with COVID-19. Taking the data deposited in the GISAID database, we made predictions using the IEDB analysis resource (TepiTool) to gauge how variants in the SARS-CoV-2 genome may change peptide binding to the most frequent MHC-class I and -II alleles in Africa, Asia and Europe. We caracterized how a single mutation in the wildtype sequence of of SARS-CoV-2 could influence the peptide binding of SARS-CoV-2 variants to MHC class II, but not to MHC class I alleles. Assuming the ORF8 (L84S) mutation is biologically significant, selective pressure from MHC class II alleles may select for viral varients and subsequently shape the quality and quantity of cellular immune responses aginast SARS-CoV-2. MHC 4-digit typing along with viral sequence analysis should be considered in studies examining clinical outcomes in patients with COVID-19.